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Abstract: We quantify uncertainties in the Monte-Carlo simulation of inclusive and dijet 
final states, which arise from using the MC@NLO technique for matching next-to- 
leading order parton level calculations and parton showers. We analyse a large 
variety of data from early measurements at the LHC. In regions of phase space 
where Sudakov logarithms dominate over high-energy effects, we observe that the 
main uncertainty can be ascribed to the free parameters of the parton shower. 
In complementary regions, the main uncertainty stems from the considerable 
freedom in the simulation of underlying events. 



1 Introduction 



QCD jet production constitutes an important background in a variety of searches for theories beyond the 
Standard Model [l}{6]. At the same time, measurements of inclusive jet and dijet cross sections are used 
to constrain parton distributions [7l|9] and to determine the value of the strong coupling 10 12 . Despite 



the tremendous importance of QCD jet production, precise predictions of event rates and kinematics using 
higher-order perturbation theory remain challenging. Only up to four-jet final states have been computed 
at the next-to-leading order so far 13-20 . Phenomenologists therefore typically rely on the simulation of 



high-multiplicity signatures by Monte-Carlo event generators. 

The ATLAS and CMS experiments at the CERN Large Hadron Collider have recently measured inclusive jet 



and dijet production 21 - 27 , with many observables implicitly probing higher-order effects. The outstanding 



quality of these data allows to validate and refine existing Monte-Carlo tools. The scope of this publication 
is a quantification of related perturbative and non-perturbative uncertainties. 

Calculating next-to-leading order QCD corrections to arbitrary processes has become a highly automated 
procedure, limited only by the capacity of contemporary computing resources. Infrared subtraction tech- 
niques 28 -30 are implemented by several general-purpose matrix element generators 31-35. The compu- 



tation of virtual corrections is tackled by a variety of dedicated programs 36 -51 . Turning the parton level 



result into a prediction at the particle level then requires a matching to the parton shower in order to imple- 
ment resummation. Two methods have been devised to perform this matching procedure, the Mc@NLO 52 
and the POWHEG 53 54 technique. While both are formally correct at the next-to-leading order, they 



exhibit subtle differences, which have been in the focus of interest recently 55 56 . We shall continue this 



study to some extent and perform a detailed comparison of scale uncertainties with resummation uncertain- 
ties as well as ambiguities arising from the Monte-Carlo simulation of the underlying event. We will employ 
the Mc@NLO technique to match next-to-leading order parton-level results for dijet production with the 



parton shower as implemented in the event generator SHERPA 57 58 . Virtual corrections are obtained from 
the BLACKHAT library 20 41 59 . Earlier studies of inclusive jet and dijet production used the POWHEG 



approach 60 . They exhibit a large dependence on the parton-shower and underlying-event model |21|. We 



expect that the conclusions drawn from our study will also apply to the simulation provided in [60] , as the 
MC@NLO and POWHEG techniques are of the same formal accuracy. 

The outline of this paper is as follows: Section[2]introduces the theoretical framework for our study, including 
a description of the new developments in SHERPA, which allow to perform the variation of scales in a manner 
consistent with analytical resummation techniques. Section [3] presents results and discusses the size and 
relative importance of the various types of uncertainties. Section [4] contains some concluding remarks. 



2 The Mc@Nlo matching method and its uncertainties 



This section outlines the essence of the Mc@NLO technique for matching next-to-leading order matrix el- 



ements and parton showers. We follow the notation introduced in 55 and report on an extension of the 



MC@NLO implementation therein, which allows to vary resummation scales in a more meaningful way. We 
point out the free parameters of the Mc@NLO method, which will be used to obtain quantitative predictions 
in Sec. [3l 



Notation 

In the following, B($b) will be used to label Born squared matrix elements, defined on the Born phase space 
which are summed/averaged over final-/initial-state spins and colours and include parton luminosities 
as well as symmetry and flux factors. Squared matrix elements of real emission corrections are denoted 
by R($a). They are defined on the real-emission phase space $r. Virtual corrections, including collinear 
countcrtcrms, are denoted by V. Real and virtual corrections induce infrared singularities of opposite sign, 
which cancel upon integration 61 62 . In order to exploit this cancellation for the construction of Monte- 



Carlo event generators, subtraction formalisms are invoked 28 29 , which introduce real subtraction terms 

(S) . . . (s) — — 

DJ - 1 and their corresponding integrated counterparts I ' 



ij,k — --t— ™— o ™,~ ~— - „„_~ if -™„„ 

In the MC@NlO method, one defines additional Monte-Carlo countcrtcrms 



D 



(A) 



y k , which represent the evo- 
lution kernels of the resummation procedure. These counterterms must necessarily have the correct infrared 
limit in order for the method to maintain full NLO accuracy 55 . Therefore, full colour and spin information 



needs to be retained. Ordinary parton-shower evolution kernels are recovered from the MC counterterms by 
taking the limit N c — > oo and averaging over spins. 

Away from the collinear limit, the form of MC counterterms is less constrained, which essentially presents a 
source of uncertainty of the Mc@NLO method. This particular uncertainty will not be addressed here as it 



can be reduced by matrix-element parton-shower merging at the next-to-leading order 63 - 65 



It is particularly useful to identify the MC counterterms with infrared subtraction terms, D 



(A) 
ij,k 



D 



(S) 



up 



to phase space constraints. This procedure was advocated in [55] , but the corresponding implementation 
suffered from unknown integrals in the integrated subtraction terms in the case that the phase-space was 
parametrised in terms of parton-shower evolution and splitting variable. In this publication, the problem is 
solved by performing the integral of the remainder term D 



(A) 
ij,k 



Djj numerically. 
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The McONLO method 



In terms of the above defined quantities, omitting flavour- and phase-space mappings, the expectation value 



of an arbitrary infrared safe observable O in the Mc@NLO method is given by 52 
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Therein, Born phase space configurations <&b are assigned next-to-leading order weights according to 



B^(<S> B ) = B($ B )+V($ B ) + £ i^(*b) 
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(2.2) 



The real-emission phase space associated with parton emission off an external leg tj of the Born configuration 
$b can be factorised as <I>_r = $b • I n this context, k denotes spectator partons in the splitting 

ij, k — > z, j, k, which are used to absorb the recoil when a splitting parton ij is put on-shell in the subtraction 



procedure. The emission phase space, d<&oj B , can be parametrised as d$p|p cx dtdzdcj), i.e. in terms of 



ij.k 



R\B 



R\B 



k^. In the above equations, we always assume t 



an evolution variable t, a splitting variable z and an azimuthal angle 
identified with some transverse momentum l ~ 

t($R,<S> B )- 

We call Eq. (|2.2|) the next-to-leading order weighted Born cross section. The sum over parton configurations 



The evolution variable t is usually 

L \^R\B> 



and the integral over the emission phase space in Eq. (2.2) are both evaluated using Monte-Carlo methods, 

US) 



while keeping the Born phase space point, $5, fixed. Because of the choice d| A ^ = D^ ; fe , the integrand 
varies only mildly. A single phase-space point is therefore sufficient to obtain a reliable estimate of the 

I 1 • • • (A) 

integral. Note that in 55 the corresponding integral was absent as the phase-space constraints on ■ k 
— (s) 

were chosen to be the same as on D > . ' k . While this method lead to fewer fluctuations of the MC integral, it 
severely restricted the flexibility of the Mc@NLO and hampered the correct assessment of uncertainties. In 
this publication we are able to lift this restriction while still maintaining full next-to-leading order accuracy 
of the simulation through incorporating full colour and spin information in the resummation. 

The resummation procedure itself is encoded in the square bracket multiplying the NLO-weighted Born 
matrix element, & A \ on the first and second line of Eq. (2.1). Note that the square bracket is unitary by 
construction. The overall Sudakov factor A^ A ) is defined as 



(2.3) 



with the partial Sudakov factors given by 



A$(M0 = exp<- 



: -q,g t 
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Again, t = is implied and the dependence of the Sudakov factor on the Born phase space config- 

uration is implicit. This Sudakov factor differs from the ordinary parton-shower Sudakov by including full 
colour and spin correlations. Its implementation is detailed in 55 . The factors Sb, Sr and Sij account for 



the potentially different symmetry factors present in the Born and real-emission matrix elements and the 



3 



parton-shower expression, respectively. In the latter, identical particles produced at different scales t are 
distinguishable, leading to a factorisation of symmetry factors along the evolution chain. The third line in 
Eq. (2.1) encodes the non-logarithmic remainder terms of the next-to- leading order real-emission correction. 
Subsequent parton-shower evolution is effected on both terms respecting the emission scales already present. 
If necessary, truncated parton-shower emission are inserted 
parton shower. 



53 to retain the logarithmic accuracy of the 



Uncertainties 



The evolution variable in the Monte-Carlo counterterms is limited from above by the resummation scale 
squared, /ig. This scale was introduced in the context of analytic resummation 66 -70 . It can be used to 
assess the uncertainties associated with the resummation programme. We will make extensive use of this 
possibility in Sec. [3] 

When defining scales for different parts of the calculation, there are certain restrictions that have to be 
adhered to. In principle, both the factorisation scale, fip, and the renormalisation scale, (Ir, in the NLO 
weighted Born ME and the hard remainder function can be chosen freely. The difference induced by different 
scales is formally of 0(a 2 s ) relative to the Born contribution. It can, however, have a sizable impact in 
practice. 

A different scale can be chosen for the resummation kernel in the square bracket of Eq. (2.1), which corre- 



sponds to the scale employed by the parton shower. This scale is required to be consistent with the one used 
in the shower itself. In order to achieve full next-to-logarithmic accuracy, it must be of the functional form 
of the relative transverse momentum of the splitting, kj^ , 



71 



as k_i_ — s- 

Additional uncertainties arise from subsequent parton showers, hadronization and the simulation of multiple 
parton scattering. We comment on these effects in the following section. 



3 Results 



In this section results generated with the MC@NLO algorithm detailed previously are presented for inclusive 
jet and dijet production. Monte-Carlo predictions are compared to a wide variety of measurements made 
by both the ATLAS and CMS experiments at the LHC at 7 TeV. The automated implementation of the 
MC@NLO algorithm in the Monte-Carlo event generator SHERPA, detailed in 55 , is used where the only 



non-automated ingredient, the one-loop matrix clement, is interfaced from BLACKHAT 15 16||20| employing 



the methods of |72|. Further QCD evolution is effected using SttERPA's built-in CSS parton shower 

hadronization corrections 
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Non-perturbative corrections, including multiple parton interactions 74 



75 



76 



and hadron decays 77 , are calculated using phcnomcnological models tuned to data. The standard tune for 
SHERPA-1.4.0 has been used. Soft-photon corrections are simulated using 78 . The CT10 parton distribution 



functions 79 have been employed throughout. All results are presented at the particle level, using only stable 



final state particles with a lifetime longer than 10 ps. 

For the central theoretical prediction the scales have been chosen as 



A*F = Mil = 4 Ht 



and 



MQ 



Therein, Ht is defined as the sum of the scalar transverse momenta of the jets found in the partonic process 
before applying any resummation. These partonic jets are defined using the anti-fcj^ algorithm [80,81 with 
R = 0.4 and p™" 1 = 20 GeV. p± is the transverse momentum of any of the two jets of the Born phase 
space configuration upon which the MC@NLO procedure is effected. To estimate the intrinsic uncertainty 
of the predictions the perturbative scales /if and (Ir have been varied independently by the conventional 
factor of 2 around the central scale while the resummation scale is kept fixed. For the resummation scale 
variation, taking into account the simple form of the exponent of the Sudakov factor, the prescription of 67 
is followed, varying /ig by a factor of y2 around the central choice while both fip and /jr are kept fixed. 
Further, non-perturbative uncertainties, i.e. the impact of shortcomings of the phenomenological models, 
have been assessed following the prescription outlined in 82 : Alternative tunes increasing and decreasing 



the mean charged multiplicity in the transverse region by 10% were used to estimate the uncertainty in the 
multiple parton interaction modeQ Exchanging the cluster hadronization model of 75 76 for the Lund 



1 This corresponds to changing the switch SIGMA_ND_FACT0R=f0.03 (increase/decrease) around the central tune value. 
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Figure 1: Colour scheme used to display various uncertainties. Overlapping uncertainties will be displayed 
by adding the colours as indicated. 
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Figure 2: Inclusive jet cross section compared to ATLAS data 
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83 has been found to have negligible impact on jet observables in previous 



string hadronization model 

studies [55j[82], and is therefore not considered here. 

All observables studied in the following have been calculated using the same event sample. It is defined by 
requiring at least two anti-fcj_ jets (R — 0.4) with p± > 10 GeV at the parton level before applying any 
resummation, of which at least one must have p± > 20 GeV. Fig. [T] presents the colour scheme used to 
display the individual uncertainties and their overlaps. 



Inclusive jet rates 

The first observables to study are inclusive jet production rates. These have been measured by the ATLAS 
collaboration 121]. Jets are defined at the particle level using the anti-fcj^ algorithm with R = 0.4 and p± > 60 
GeV within \y\ < 2.8. Jets are ordered in transverse momentum. Additionally, the leading jet is required to 
have p ± > 80 GeV. 

Fig. [2] presents the results. We observe good agreement between our Monte-Carlo simulations and experi- 
mental data. The renormalisation and factorisation scale uncertainty amounts to approximately 7% for the 
dijet inclusive cross-section, which is described at next-to-leading order accuracy, while it increases to 14% 
for the three-jet inclusive rate, which is described at leading order accuracy only. All higher multiplicity jet 
inclusive rates are described at the logarithmic accuracy of the parton shower only and therefore inherit the 
scale uncertainty of the inclusive three-jet rate. The resummation uncertainties, indicating the observables 
sensitivity to multiple higher-order soft emissions below the resummation scale, is slightly larger: 8% for 
the dijet inclusive cross section and 35% for the three-jet inclusive rate. They steadily increase for higher 
jet multiplicities. The non-perturbative uncertainties, on the other hand, are negligible, contributing from 
^0.2% for the dijet cross section to ~6% for the inclusive 5 jet cross section. 
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Jet transverse momenta (anti-kt R-0.4) 
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Figure 3: Jet transverse momenta compared to ATLAS data 

Scalar sum of jet transverse momenta (anti-kt R-0.4) 
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Figure 4: Scalar sum of jet transverse momenta compared to ATLAS data 21 
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3 jets over 2 jets ratio (anti-kt R-0.5) 
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Figure 5: 3-jet over 2-jet ratio in dependence on the scalar transverse momentum sum of all (CMS) and 



the two leading (ATLAS) jets in comparison to CMS [24] and ATLAS data [21 . 

Jet transverse momenta 



The same analysis 21 studied also the pj_-spectra of the individual jets. The event selection is the same 
as above, except that subleading jets with p± > 60 GeV may or may not be present in case of the leading 
jet p^. Fig. [3] displays the results compared to ATLAS data. Again, we observe good agreement with 
our Monte-Carlo predictions. The two leading jets' transverse momenta, both calculated at next-to-leading 
order accuracy over large parts of the phase space, show the characteristically small renormalisation and 
factorisation scale dependencies of ^5-10%. Their rcsummation scale dependence is comparably very small 
throughout, ranging from ^5% at low p± to ^20% at large p±. This is expected from all choices of 
the resummation scale being smaller than the second jet's transverse momentum. Any influence therefore 
stems from the mismatch of (Mc@Nlo) parton shower evolution and jet reconstruction. The transverse 
momentum of the third jet, being calculated at leading order accuracy, and the fourth jet, determined at 
leading logarithmic accuracy only, exhibit much larger scale variation, both for the renormalisation and 
factorisation scales and the resummation scale. Non-perturbative uncertainties are small in comparison. 
Similarly, Fig. [4] displays the scalar sum of the individual jet transverse momenta in events with at least 
two, three or four jets. For these observables good agreement is found as well. The perturbative and 
non-perturbative uncertainties are comparable to those of the individual jet transverse momenta. 

3-jet over 2-jet ratio 

The next observable to be examined is the relative rate of inclusive three-jet events compared to inclusive two- 
jet events. The CMS collaboration measured this ratio in dependence on the scalar sum of all jet transverse 
momenta, Ht, [24] . Within this analysis events with at least two and three jets, respectively, defined with 
the anti-fcj_ algorithm with R = 0.5, p± > 50 GeV, \y\ < 2.5, and H T > 0.2 TeV were selected. The 3-jet 
over 2-jet ratio is then defined as R32 = (d<T>3j e t/d7?T)/(d(7>2jet/di?T)- The comparison of the presented 
calculation to data is shown in the left panel of Fig. [5] Good agreement between MC predictions and data is 
observed. The scale variations in calculating (d<7>3j e t / dHx) , described at leading order, and (dcr>2jot/di?T), 
described at next-to-leading order, were done simultaneously because both observables were calculated from 
the same event sample. Consequently, the renormalisation and factorisation scale uncertainty are large at 
small Ht but largely cancel at large Ht- The resummation uncertainty behaves similarly, but shows an 
opposite asymmetry at small Ht- The non-perturbative uncertainties are much smaller for small Ht but 
grow to equivalent size in the large- Ht region. 

The analysis of 21 also studied the 3-jet over 2-jet ratio both as a function of the the scalar transverse 



(2) 

momentum sum of the two leading jets, H T , and as a function of the transverse momentum of the leading 
jet, p 1 ^, only. The results are displayed in the right panel of Fig.JHjand in Fig. [6] respectively. Both analyses 
show the same level of agreement between MC predictions and data. Scale uncertainties and non-perturbative 
uncertainties are of similar size as in the CMS analysis. 
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Figure 6: 3-jet over 2-jet ratio in dependence on the leading jet transverse momentum in comparison to 



ATLAS data 21 . 
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Figure 7: Azimuthal decorrelations of the leading jets compared to CMS data [25]. The individual p 1 ^ 
ranges on the left-hand side plot are ordered as indicated in the right-hand side ratio plots 
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Inclusive jet transverse momenta in different rapidity ranges 
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Figure 8: Inclusive jet transverse momenta spectra different rapidity ranges compared to ATLAS data 

The different rapidity ranges on the left-hand side plot are ordered as indicated in the right-hand 
side ratio plots, being multiplied, from top to bottom by factors of 1, 3-10 -2 , 10~ 3 , 3T0 -5 , 10~ 6 , 
3-10~ 8 , and 10 -9 , for readabilities sake. 



Azimuthal decorrelations 

Next, the correlations between the two leading jets are examined. The CMS collaboration measured the 
dijet azimuthal decorrelations, i.e. the A(j> separation of the two leading jets, in [25] . Therein, the jets are 
defined using the anti-fcj^ jet algorithm with R = 0.5 and a minimum transverse momentum of p± > 30 GeV 
within a rapidity interval of \y\ < 1.1. Events with at least two such jets were classified according to the 
leading jet's transverse momentum into five mutually exclusive regions: p^ ad g [80,110] GeV, [110,140] 
GeV, [140, 200] GeV, [200, 300] GeV, and [300, oo) GeV. The present calculation provides next-to-leading 
order accuracy at Atfi — ir, leading order accuracy in the region |7T < A(j> < it, and leading logarithmic 
accuracy in the region A0 < |7r. The results are shown in Fig. [TJ Good agreement between data and MC 
prediction is found. The renormalisation and factorisation scale uncertainties are accordingly small in the 
region A(f> > |7r and increase towards lower A(f> values. Varying the resummation scale leads to similarly 
large uncertainties, while non-perturbative uncertainties only play a minor role. However, by normalising 
the observable to the inclusive dijet cross section, the scale uncertainties are artificially reduced. 



Inclusive jet transverse momenta in different rapidity ranges 

To further study correlations between multiple jets produced, it is useful to consider double-differential 
observables studied by the ATLAS collaboration [22]. We start with the inclusive jet transverse momentum 
in different rapidity ranges. Jets are defined using the anti-fcj_ jet algorithm with R = 0.4, p± > 20 GeV 
and \y\ < 4.4. Every jet is considered in the analysis. The contribution from the first two jets in the region 
where at least two jets are present is described at next-to-leading order accuracy, while the contribution 
of a possible third jet is described at leading order. All contributions of subsequent jets are described at 
leading logarithmic accuracy only. Thus, the overall accuracy of these observables is a mixture of the above. 
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Dijet invariant mass spectra in different rapidity ranges 
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Figure 9: Dijet invariant masses in different rapidity ranges compared to ATLAS data 22 . The rapidity 
different difference ranges on the left-hand side plot are ordered as indicated in the right-hand 
side ratio plots, being multiplied, from top to bottom by factors of 10 12 , 3 TO 10 , 10 9 , 3 TO 8 , 10 6 , 
3T0 4 , 10 3 , 30, and 1, for readabilities sake. 



Fig. [8] shows the result of the presented calculation compared to data. The agreement in all but the most 
forward rapidity ranges is good. Renormalisation and factorisation scale uncertainties are small for central 
jet production while they grow larger with increasing rapidity. The resummation scale uncertainty behaves 
similarly albeit being larger in magnitude throughout. Non-perturbative uncertainties are small, except in 
the very forward region, close to the beams. At very large rapidities the transverse momentum is no longer 
a good measure of the hardness of the process 13 14 84 . Instead a scale taking into account the dijet 



invariant mass should be used. Such consideration applied to an i/^-based scale, taking into account real 
emission dynamics, is proposed to take the following form 

of \y-yhao B t\ 



Vr/f 



i H (y) 

4 n T 



E 

iGjets 



\p 



(3.1) 



wherein yboost = l/ n jct ■ Siejets ^ ne rapidity of the n-jet system. The factor / is chosen to suitably 
interpolate between the invariant-mass-like behaviour and transverse- mome ntum- like behaviour. For / = 0.3 

o 0.3y* ~ 



13U14 
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and the presence of only two jets it reduces exactly to the scale proposed in 
"112/(4 cosh(0.7y*)). This scale choice, however, is only beneficial for describing Hata in this and the following 
analysis and either shows no impact or even reduces the agreement with the observed experimental data in 
all other analyses considered in this publication^] It is therefore not adopted as a central scale choice. 



Dijet invariant masses in different rapidity ranges 

Another doubly differential observable studied in [22] is the dijet invariant mass. Events with at least two 
anti-fc_L jets with R = 0.4 and p 1 ^ > 30 GeV and p ™ blcad > 20 GeV within \y\ < 4.4 are considered. 

2 Taking /j.r/p = j as the central scale leads to an underestimation of the inclusive three-jet rate presented in Fig. [2] 

by 30%, for example. 
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The dijet invariant mass is defined m\i = \J (p lcad + j,subicad)2 anc j bi nnec j j n to various mutually exclusive 
ranges of y* — | |yi oa d — ysubleadli the rapidity separation of the two jets. Fig. [9] displays MC results 
compared to ATLAS data. For small rapidity separations the agreement with data is good. At large y* the 
renormalisation and factorisation scale is too low, resulting in too large MC predictions. This situation is 



improved by choosing the scale Eq. (3.1 1 instead. The renormalisation and resummation scale uncertainties 
are small at small y* and increase towards larger y* . They also remain approximately constant over the whole 
considered to 12 range. Resummation scale uncertainties, on the other hand, are larger and display a definite 
TO12 dependence. This can be seen as an indication that high-energy resummation, which goes beyond the 



collincar limit used in the parton shower, becomes important 85 - 87 . Non-perturbative uncertainties are 
only non-negligible at low invariant masses or large rapidity separations, when at least one jet is likely to be 
close to either of the two beams. 

As discussed in the previous paragraph, taking \ Hjf' as the central scale improves the description of this 
observable but leads to a worse description of the other observables investigated in this paper. It is therefore 
not adopted as the central scale. 

Gap fractions 



A different way to probe the radiation pattern was explored by the ATLAS collaboration in 23 . Therein, 
events were selected containing at least two jets, defined using the anti-/c_L algorithm with R = 0.6, each with 
p± > 20 GeV within y < 4.4. Within these events a dijet system is then identified using either the two largest 
transverse momentum jets (leading jet selection) or the widest separated jets (forward backward selection). 
For both definitions an average transverse momentum p± = ^(p^ 1 +P^ t2 ) of at least 50 GeV is required. 
To characterise the subsequent radiation pattern two variables are used. The gap fraction, i.e. the fraction 
of events that do not exhibit any further radiation above some Qq within the Ay rapidity range spanned 
by the dijet system, and the mean number of jets with p± > Qq m the same Ay region. Fig. [10] displays 
a comparison of MC results with data for the gap fraction (Qq = 20 GeV) in dependence on the rapidity 
separation of the dijet system and its p± for both selections. The agreement is good throughout the probed 
region. In both selections the renormalisation and factorisation scale uncertainty is the leading uncertainty 
at low p± while the resummation uncertainty dominates at large p±. Non-perturbative uncertainties are 
generally larger than for most other observables considered in this paper, increasing with Ay and p±. In 
case of the forward backward selection they are of the same magnitude as both perturbative uncertainties 
in the large Ay and p± region. 



Similarly, Fig. 11 displays the average number of jets (Qq = 20 GeV) in dependence on the rapidity separation 
of the dijet system and its p± for both selections. Again, good agreement is found. All uncertainties are 
small and of comparable size for small Ay throughout the p± range, and steadily increasing for larger Ay. 
While a resummation scale variation produces largely p± independent uncertainties the renormalisation and 
factorisation scale uncertainties are larger for small average transverse momenta than for large ones. The 
non-perturbative uncertainties show the opposite behaviour. 

Event shapes 

Traditional observables not being described without resummation are event shapes. The CMS collaboration 
measured the central transverse thrust and the central transverse thrust minor in multijet production in [26| . 
The sample is defined by requiring at least two jets, defined using the anti-fcj_ algorithm with R = 0.5, 
p± > 30 GeV and |f?| < 1.3. The selected events are then categorised into three mutually exclusive regions 
according to the leading jet transverse momentum. The observables are defined as follows 

EiejetelP-M-^H , m Eiejots IP-L.i X h T,C I 

t±,c = 1 — max — — — ; — and J m ,c — 



^ EiGjets "' E lG jcts W-L,i\ 

Therein, the vector n*r,c is defined as the vector minimising t_l,c. Only jet momenta are taken into account. 



The results of the presented calculation are compared to the experimental data in Fig. 12 Good agreement 
between MC predictions and data is found. The dependence of both observables on the renormalisation and 
factorisation scale, despite being calculated at most at leading order accuracy, largely cancels due to their 
normalisation. However, they do show a large dependence on the resummation scale, as expected. This 
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Figure 11: Mean jet number in gap in dependence of mean transverse momentum and rapidity separation 



of dijet pair for both selections compared to ATLAS data 23 . 
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Central transverse thrust in different leading jet p± ranges 
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Figure 12: Central transverse thrust and central transverse thrust minor compared to CMS data 26 . 
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Forward energy flow in dijet events, p ± > 20 GeV 
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Figure 13: Forward energy flow compared to CMS data 27 



is largest when rj_,c C?m,c) is close to zero (one). The non-perturbative uncertainties have the opposite 
behaviour. 



Forward energy flow 



The last observable to be studied is the forward energy flow as measured by the CMS collaboration 27 . In 
this analysis, events with at least two jets, defined by the anti-fc_L algorithm with R — 0.5 and p± > 20 GeV 
are required. The two leading jets are further required to lie within \rj\ < 2.5 and satisfy | A<fi— ir\ < 1, i.e. to 
produce a nearly back-to-back topology. Within this event sample the energy flow, defined as average energy 
E per event per pseudo-rapidity interval dry, is calculated and compared to the measured data. Despite 
a small difference in shape good agreement is found. Again, due to the normalisation of the observable, 
renormalisation and factorisation scale uncertainties as well as resummation scale uncertainties are small. A 
comparably large uncertainty stems from the non-perturbative modelling uncertainties, ranging up to ^10%. 
This is not unexpected since in this very forward region, close to the beams, non-factorizable components of 
the inclusive cross section play a large role. 



4 Conclusions 

We have presented a detailed analysis of uncertainties associated with the simulation of inclusive jet and dijet 
production using methods for matching next-to-leading order QCD calculations and parton showers. We 
have analysed factorisation and renormalisation scale dependence as well as variations originating from the 
choice of resummation scale. We have compared to uncertainties originating from the freedom in choosing 
parameters in the Monte-Carlo simulation of multiple parton scattering. 

These three types of uncertainties represent different degrees of freedom in the Monte-Carlo simulation: While 
the renormalisation and factorisation scale dependence probe the impact of higher-order QCD corrections 
to the hard process, the resummation scale dependence quantifies, to some extent, uncertainties related to 
parton evolution. Variations of the MPI tune are used to estimate uncertainties related to non-perturbative 
dynamics. 

The results of our Monte-Carlo simulation have been compared to a variety of data taken by the ATLAS 
and CMS experiments at the CERN LHC. Good agreement is found for almost all observables. Exceptions 
are the inclusive jet transverse momenta at large jet rapidity and the dijet invariant masses at large average 
rapidity. Discrepancies are attributed to the choice of renormalisation and factorisation scale, which is given 
by one quarter of the visible transverse energy. Despite taking into account real-emission dynamics, this 
scale does not give a realistic measure of the hardness of events at large individual jet rapidities. A modified 
scale, with jet transverse momenta weighted by their rapidity w.r.t. the centre of the partonic system, leads 
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to better agreement in the forward region, but deviations are observed for central jet production. 

Uncertainties related to higher-order corrections to the hard process are most important for exclusive multi- 
jet final states. Similarly, uncertainties related to the resummation procedure are most significant in the 
region where jet production is modelled by the parton shower. We expect that both uncertainties can be re- 



duced by application of matrix-element parton-shower merging methods at the next-to-leading order 63 - 65 



A corresponding analysis is forthcoming. A different role is played by non-perturbative uncertainties. They 
are most significant in regions where (semi-)soft particle production dominates over multi-jet effects. They 
can be reduced only by better constraints on the non-perturbative dynamics through additional measure- 
ments of underlying event activity and particle flow. 

Our analysis should help to better understand the quality of Monte-Carlo predictions, which are obtained 
using matching methods like Mc@NLO and POWHEG. Ascribing reliable uncertainties to such predictions 
will remain an important task in the immediate future of LHC physics. 
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